Optical Lattice Induced Light Shifts in an Yb Atomic Clock 
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We present an experimental study of the lattice induced light shifts on the So 
clock transition {fciock ~ 518 THz) in neutral ytterbium. The "magic" frequency, u mag i C , 
174 Yb isotope was determined to be 394 799 475(35)MHz, which leads to a first order light shift 
uncertainty of 0.38 Hz on the 518 THz clock transition. Also investigated were the hyperpolarizability 
shifts due to the nearby 6s6p 3 Po — > 6s8p 3 Po, 6s8p 3 P2, and 6s5/ 3 f2 two-photon resonances at 
759.708 nm, 754.23 nm, and 764.95 nm respectively. By tuning the lattice frequency over the two- 
photon resonances and measuring the corresponding clock transition shifts, the hyperpolarizability 
shift was estimated to be 170(33) mHz for a linear polarized, 50 /iK deep, lattice at the magic 
wavelength. In addition, we have confirmed that a circularly polarized lattice eliminates the J = 
— > J = two-photon resonance. These results indicate that the differential polarizability and 
hyperpolarizability frequency shift uncertainties in a Yb lattice clock could be held to well below 

io- 17 . 



PACS numbers: 32.10.Dk, 06. 30. Ft, 32.70.Jz, 39.30. 



The latest generation of atomic clocks are based on 
narrow optical transitions in atoms or ions. Moving from 
microwaves to optical frequencies greatly increases line 
quality factors, and is enabled by methods of laser cool- 
ing and trapping, highly stabilized CW laser sources, and 
the rapid advances in optical frequency comb technology 
that allow high accuracy comparison of frequencies rang- 
ing from the RF to the UV. Optical lattice clocks, by 
confining atoms to the Lamb-Dicke regime, have the po- 
tential to achieve very high accuracy (systematic uncer- 
tainty less than 10~ 17 ), while at the same time achieving 
high stability through the use of large numbers of atoms. 

Although most initial lattice clock studies have focused 
on Sr P, 3, [H, HI j the alkaline earth- like Yb is also an ex- 
cellent candidate for an optical lattice clock 0, d, 0| • In 
addition to being relatively simple to laser cool and trap, 
Yb has a spin- 1/2 isotope ( 171 < 173 Yb, I = 1/2,5/2 re- 
spectively) and works well for magnetically-induced spec- 
troscopy of the spin zero isotopes [H, either of which can 
be used to reduce experimental complexities due to the 
nuclear substructure. In addition to clocks, Yb trapped 
in optical lattices is of interest for studies of Boson- 
Fermion interactions and quantum information (Toj . 
In this Letter, we report the first accurate measurement 
of the magic wavelength for the 1 Sq — » 3 Pq clock tran- 
sition of Yb. We also find that the hyperpolarizability 
shifts caused by two-photon transitions nearly resonant 
with the magic wavelength are comparable to those for 
Sr when scaled in terms of trap depth [i"ll|. These mea- 
surements are essential to building an accurate optical 
lattice clock. The lattice induced light shifts of our Yb 



clock can currently be determined to < 8 x 1CP 16 in frac- 
tional frequency, and there exists no practical limit to 
reducing this to the 10~ 18 level. 

A major concern for the accuracy of optical lattice 
clocks is the ability to cancel the large light shifts created 
by the trapping potential of the lattice. Practical optical 
lattice depths for Yb are 50 to 500 E r (E r is the lattice 
recoil energy), which means the corresponding light shifts 
of 100 kHz to 1MHz must be canceled to a high degree. 
The light shift, Aw, on the clock transition due to the 
lattice can be described by a power series expansion in 
the electric field strength E: 



Aw = a (wj,e) E 2 + 7 (oj h e) E 4 



(1) 



where a and 7 are the differential polarizability and hy- 
perpolarizability between the ground and excited states 
respectively, and depend on the lattice frequency, uii, and 
polarization vector e. The first term dominates, and the 
design of optical lattice clocks is such to minimize this 
polarizability term be the choice of polarization insensi- 
tive J = — > J = clock transitions and determining 
the lattice wavelength where a vanishes. At this magic 
wavelength the frequency uncertainty of the first order 
lattice induced light shift could be held below 10~ 17 

Higher order light shifts have also been shown theoret- 
ically to be below the 10~ 17 level, barring near resonant 
terms in the hyperpolarizability from nearby two-photon 
transitions However, in Yb the magic wavelength is 
in the vicinity of three two-photon resonances from the 
6s6p 3 Po excited clock state to higher lying 6s8p 3 Pj and 
6s5/ 3 i<2 states, with the closest a third of a nanometer 
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FIG. 1: Relevant energy levels of Yb, including expanded view 
of states nearly resonant through two-photon excitation at the 
magic wavelength. The two-photon detunings are expressed 
in terms of the one photon frequency. 



away (see Fig. [T]). Our simple estimate of the strength of 
this 3 Po 3 -Po transition suggested that the hyperpo- 
larizability shift could have been as large as one hertz at 
the magic wavelength. Until now, detailed studies of lat- 
tice induced light shifts in Yb have not been performed, 
and the magic wavelength had been only roughly located 
to near 759.35 nm 

Most of our current work utilizes the 174 Yb isotope and 
the magnetically induced spectroscopic (MIS) method 
■ After laser cooling and trapping on both the strong 
1 Pi singlet line and the weaker 3 Pi intercombination line, 
we load roughly 10 000 atoms into an optical lattice at 
temperature of » 15 /iK. The ID lattice is formed in the 
horizontal plane by focusing approximately 1 W of light 
to a waist of rs 30 /im, then retroreflecting the beam with 
a curved mirror. This produces a standing wave lattice 
consisting of a series of pancake shaped wells with a trap 
depth of ~ 50 yK (500 E r ). After the preparation stage, 
a static magnetic field (~ 10 mT) is turned on to mix the 
upper 3 Po clock state with the 3 Pi state, which provides 
a nonzero excitation probability between the two states 
(MIS method). While the bias field is on, the atoms are 
probed along the lattice direction with a frequency stabi- 
lized laser at the clock wavelength of 578.42 nm. Excita- 
tions to the upper clock state are detected by measuring 
the depletion of the ground state population after the ex- 
citation pulse. With a field of s=s 1.5 mT and about 12 /iW 
of probe light focused to a 35 /im spot, we generate a 7r 
pulse in 64 ms. For operation as a frequency reference, 
the stable clock laser is locked to the atomic resonance 
by use of an acousto-optic modulator to remove the slow 
drift of the reference cavity. Further details of the exper- 
imental apparatus can be found in reference [12[ ■ 

For the measurements of the magic wavelength, we 
compared our clock frequency shifts to the Ca clock [l3j| 
using an octave spanning Ti:Sapphire optical frequency 
comb [14( | . The Ca clock runs in a simplified mode that 
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FIG. 2: a) The shift of the Yb optical clock frequency versus 
lattice power for twelve separate lattice wavelengths. As the 
lattice wavelength is increased, the differential polarizability 
(slope of the lines) goes from large and negative, through 
zero at the magic frequency, v m agic = 394 799 475(35) MHz, 
(A m(l£ ,i C » 759.3537 nm), to large and positive. The data are 
individually offset so the fit lines pass through the origin, 
b) Differential polarizability plotted as a function of lattice 
frequency. The fit line was used to interpolate the magic 
frequency with an uncertainty of 35 MHz. 



provides high stability (averaging down to 3 x 10 -16 frac- 
tional frequency instability in w 200 s) and high reliabil- 
ity. In addition, comb was used as a reference to stabilize 
the frequency of the lattice laser. 

A rough determination of the magic wavelength was 
made by observing the width and asymmetry of the 
spectroscopic lines of the clock transition. Then the 
clock frequency shift versus lattice intensity for twelve 
lattice wavelengths was measured by counting the beat 
frequency between our Yb locked laser and the comb, 
which was stabilized to the Ca clock (see FigJ5^). The 
slope of each wavelength set, which is proportional to 
the differential polarizability at that wavelength, was ex- 
tracted with a linear weighted least squares fit. These 
slopes were then plotted as a function of the lattice fre- 
quency (see Fig. Wp)- The magic frequency is v m agic = 
394 799 475(35) MHz, or X maglc « 759.3537 nm, as deter- 
mined by the zero crossing. Combined with the slope 
{da/dv m agic = — 22(1) mHz/(£y-GHz)), the estimated 
error gives an uncertainty in the shift of 0.38 Hz, for a 
lattice depth of 500 E r . This uncertainty is limited only 
by the data presented here. With additional measure- 
ments, the magic wavelength and first order light shift 
uncertainty will be continually refined and reduced by 
inter-clock comparisons. Importantly, a 100 kHz mea- 
surement uncertainty in the magic wavelength is feasi- 
ble, which would produce an clock shift uncertainty of 
less than 10~ 18 . 
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2y AC Stark Shift at 6s6p P ->6s8p P resonance 



3 P„ - 3 F„ Two Photon Resonances in 171 Yb 




FIG. 3: AC Stark shift on the clock transition with the lat- 
tice tuned near the 6s6p 3 Po *-> 6s8p 3 Po two-photon resonance 
with linear (solid points) and circular (open squares) lattice 
polarization. Extrapolating the dispersion shaped fit to the 
magic wavelength gives a shift of 200(30) mHz for a 500 E r 
trap depth. The resonance disappears in the circular polar- 
ized case due to selection rules. 
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FIG. 4: AC stark shifts on the Yb optical clock frequency 
with the lattice laser tuned near the 6s6p 3 P -> 6s5/ 3 P 2 two- 
photon resonances. The two resonances are the F — 5/2, 3/2 
hyperfine components of the 3 p2 state in the 171 Yb isotope, 
and are split by 4240(20) MHz. The fit gives a total clock 
shift of 3.5 mHz when extrapolated to the magic wavelength. 



With the first term in Eq. Q] minimized by operating 
at the magic wavelength, it is important to determine 
the size of the second term. The largest contributors to 
the hyperpolarizability are the near resonant two-photon 
transitions from the upper clock state to the 6s8p 3 Po.2 
states at 759.7082 nm (J = 0) and 754.23 nm (J = 2) 
and to the 6s5/ 3 i<2 state at 764.95 nm. (The two-photon 
transition to the 6s8p 3 P\ state at 758.41 nm is not al- 
lowed, due to selection rules.) The closest two-photon 
transition is the 3 Pq — > 3 Pn transition, which is -184 GHz 
(-368 GHz two-photon detuning) from the magic wave- 
length (see Fig. [I]). 

To measure the effect of the two-photon transition on 
the clock frequency, we tuned the lattice close to this 
resonance and measured the 1 So — > 3 Prj clock shift, sim- 
ilar to the approach used for Sr 11[. With the lattice 
detuned from the magic wavelength, the spectroscopic 
lines observed on the clock transition were slightly asym- 
metric and broadened to about 1.5 kHz. For these broad 
lines our reference cavity was sufficiently stable (400 Hz 
linear drift over a few hours), so that the Ca clock was 
not needed as a reference. 

The light shift on the clock transition with the lattice 
frequency tuned close to the 6s6p 3 Pn — ► 6s8p 3 Pn two- 
photon resonance is shown in Fig. [3J At each point of 
lattice laser frequency, we scanned over the clock reso- 
nance 6 to 8 times. The line shapes were then fit with a 
Lorentzian line shape to determine the line center. The 
measured clock frequency shifts versus lattice wavelength 
are then fit with a dispersion line shape, which when 
extrapolated to the measured magic wavelength gives a 
clock shift of 0.80(12) fiRz/E^. For our lattice depth of 



ps 500 E r , this leads to a clock shift of 200(30) mHz or 
4 x 10~ 16 fractionally. It must be noted that the asso- 
ciated uncertainty is mainly due to our determination of 
the peak lattice intensity, as the extrapolated fit supports 
a much smaller uncertainty. 

We also performed a similar set of light shift measure- 
ments with circularly polarized light. In this case, the 
3 Pn — > 3 Pn two-photon transition vanishes because of 
two-photon selection rules. Although the circular po- 
larization was imperfect (~ 10 % linear component from 
birefringence of the vacuum windows) , the resonance was 
not discernible in either the width or centers of the lines 
(see Fig. [3]). With the use of an even isotope, the first 
order light shift is very insensitive to the polarization, 
making the lattice ellipticity a free parameter that could 
be used to minimize or cancel the hyperpolarizability 

shift HE El- 

In comparison to the 3 Pn two-photon resonance the 
3 P2 and 3 P2 two-photon resonances are much further de- 
tuned, which would suggest their associated hyperpolar- 
izability shifts should be significantly smaller. To verify 
this, we switched to spectroscopy of the 171 Yb isotope 
because large light shift broadening due with the lattice 
tuned to 754.23 nm and 764.95 nm made efficient exci- 
tation of the clock transition with the MIS technique 
difficult. For the 171 Yb isotope, much larger Rabi fre- 
quencies could be generated resulting in higher contrast 
spectroscopic lines. 

The resulting light shift measurements of the clock res- 
onance near the 3 Pb two-photon resonance are shown in 
Fig. [U The separate two-photon resonances are the two 
hyperfine components caused by the 1=1/2 nuclear spin 
of 171 Yb, which are split by 4240(20) MHz (2120(10) MHz 
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for the lattice frequency) . The large error bars of the indi- 
vidual points in Fig. [4] are caused by the large linewidths 
and asymmetries of the observed lines. By fitting the 
clock shift versus lattice frequency for the two two-photon 
resonances and extrapolating to the magic wavelength, 
the total hyperpolarizability shift from the 3 F2 state is 
3.5(3.5) mHz for typical trap depths. However, because 
the fits are dominated by only a few points, we place a 
conservative upper limit of 10 mHz on the total shift for 
our 500 E r deep lattice. 

The two hyperfine components of the 3 Jo — > 3 P2 two- 
photon resonance were located at lattice frequencies of 
380480 706 MHz and 380 480602 MHz. Clock shift mea- 
surements of these resonances (not shown) give an asso- 
ciated hyperpolarizability shift at the magic wavelength 
of -38(8) mHz for a 500 E r trap. This is consistent with 
a calculation of the shift using knowledge of the 3 Pq res- 
onance shifts and angular momentum considerations. 

Totaling 170(33) mHz for a 500 E r lattice, the com- 
bined magnitude of hyperpolarizability shifts is non- 
negligible but manageable. Indeed, the size of the shifts 
are comparable to those of Sr when scaled in terms of E r . 
The shifts are more significant in our case than for the 
present Sr lattice clocks because we perform spectroscopy 
with considerably higher lattice intensities. However, we 
anticipate using significantly shallower traps in the fu- 
ture, which with the intensity squared dependence of the 
hyperpolarizability shift will greatly reduce the magni- 
tude of the shifts. Lattice depths of 50 E r , or w 5/uK, 
should be feasible for Yb assuming that laser cooling near 
the Doppler limit (~ 4/iK) is achieved. In addition, the 
use of vertically oriented or 3D lattices could be of ben- 
efit at lower trap depths. At 50 E r the total hyperpolar- 
izability shift would be less than 5 x 10 -18 in fractional 
frequency and with uncertainty less than 10~ 18 . Thus we 
do not anticipate the need for more precise measurements 
of the hyperpolarizability coefficients unless we choose to 
keep the trap depths above 200 E r . These hyperpolariz- 
ability shift estimates should apply to any of the isotopes 
of Yb, as the two-photon detunings are much larger than 
the isotope shifts or the hyperfine splittings. Finally, 
the relative size of these two-photon shifts allows for a 
magic ellipticity where the differential hyperpolarizabil- 
ity coefficient, 7, is zero[15|; however, good control of the 
polarization would be required. 

In this Letter, we report the reduction of the total 
lattice-induced shift uncertainty to a fractional frequency 
of less than 8 x 10~ 16 . The magic wavelength has been 
determined with an uncertainty of 35 MHz, which for our 
largest trap depths (500 E r ) produces an uncertainty of 
7 x 10 -16 fractional frequency, and further measurements 
will reduce this uncertainty. The outstanding issue of 
the possibly large hyperpolarizability due to nearby two- 
photon resonances has also been resolved. The hyper- 
polarizability shift uncertainty has been measured to less 
than 7 x 10 -17 of the clock frequency. We see no obstacle 



to reducing the lattice shift uncertainties of an Yb optical 
lattice clock to well below 10 -17 . 
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